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SPECIFICATION 
TITLE 

"METHOD, ARRANGEMENT AND COMPUTER PROGRAM PRODUCT FOR 
ANALYZING NEURONAL ACTIVITIES" 

5 BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention relates to analysis of neuronal activities in neuronal 
areas, for example of nerve structures in areas of the brain of a patient. 

Description of the Prior Art 

10 Knowledge about a mode of operation of a neuronal area and about an 

interaction of neuronal areas are fundamental to functional magnetic resonance 
tomography (fMRI), as described in A. W. Toga and J. C. Maziotta (Ed.), "Brain 
Mapping: The Methods", Ch 9: M. S. Cohen: "Rapid MRI and Functional Applications", 
Academic Press 1996, which is a further development of magnetic resonance 

15 tomography. 

Magnetic resonance tomography is an imaging method that produces 
sectional images of the human body without using damaging X-rays. 

Instead MR takes advantage of the behavior of body tissue in a strong 
magnetic field. Pathological changes in body tissue, for example in the brain or 
20 spinal cord, can be detected by this modality. 

Functional disorders in body tissue, particularly in the brain of a patient, 
however, cannot be detected by means of conventional magnetic resonance 
tomography. 

This is performed by functional magnetic resonance tomography or fMRI 
25 technology. 



Neuronal activity in areas of the brain of a patient can be measured 
indirectly by means of the fMRI technique. The BOLD (Blood Oxygenation Level 
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Dependent) signal, as it is called, is measured in individual areas of the brain, this 
signal relating to neuronal activity in the respective areas. 

Dependencies, which stem among other things from structures in the brain, 
i.e. from neuronal links between nerve cells or nerve structures, exist between the 
5 neuronal activities in the areas. 

The outcome of the fMRI measurements shows the course of activity of 
individual areas over a certain period, for example during cognitive sequences as 
a result of specific perception processes or motor tasks. 

Functional disorders, in this case in the brain, are thus inherently contained 
10 in the fMRI signals measured. 

Efficient methods for the analysis and evaluation of such fMRI 
measurements are desirable in order to be able to furnish evidence of possibly 
existing functional disorders in specific areas. 

Known methods, such as described in A.R. Mcintosh et aL, Structural 
Equation Modeling and Its Application to Network Analysis in Functional Brain 
Imaging, Human Brain Mapping, 2:2-22, 1994, are restricted to detection of 
functional relationships between various areas of the brain in certain 
predetermined tasks such as the aforementioned perception processes or motor 
tasks (functional connectivity). These functional relationships are also designated 
functional connectivity. 

In contrast to functional connectivity, however, the determination of a true 
physical connectivity, i.e. the determination of actually existing linking structures 
(of areas of the brain) independently of specific predetermined tasks, is not 
possible with these known methods. 

25 A further known method of analysis for detecting functional connectivity is 

described below. 
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The object of this known method of analysis described below is the above- 
described detection of functional relationships between various areas of the brain 
in specific perception processes or motor tasks. 

This known method of analysis is based on a predefined model of a brain, 
5 i.e. a predefined brain architecture. 

This brain architecture, predetermined from prior knowledge, defines 
general functional and/or spatial dependencies between specific areas of the brain 
in the form of a coupling matrix S, as it is called. 

The coupling matrix S has a form or structure that is fixed in accordance 
10 with the predetermined brain architecture (columns/rows) and is accordingly 
populated in certain, but not at all (matrix) positions, with changeable coupling 
strengths Si. These are changeable and are matched as part of the method of 
analysis. 

The unpopulated (matrix) positions are populated with fixed, unchangeable 
1 5 values, namely zero. 

The coupling strengths Si describe functional dependencies respectively 
between two areas of the brain or the BOLD signals measured there and 
represent the neuronal activities there. 

In this known method of analysis, the (changeable) coupling strengths Si 
20 are defined such that statistical characteristic quantities that are determined by 
this method of analysis from the fMRI measurements, can best be explained. 
Expressed differently, the probability for an occurrence of the measured data, i.e. 
the fMRI measurement or the BOLD signals, is maximized by means of the 
desired coupling strengths Si. 

25 In this method of analysis a data point s=s t represents a totality of all BOLD 

signals s1, sN of the individual n areas at a point in time t or averaged over a 
time interval t(t=[1;T]). 
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The fMRI measurement includes a large number of such data points for 
possibly differing perception processes and/or motor tasks, for which the 
corresponding BOLD signals were measured. 

In this known method of analysis, it is not the individual data points s1, s2, 
5 St, which are evaluated directly, but statistical characteristic quantities that 

emerge from these. 

For a statistical distribution of the data points s1, s2, sT it is assumed 
that it is described fully by a multivariant normal distribution, i.e. a statistical 
distribution of the first order, with a mean value m and a covariance Z: 

10 P(s\m,Z) = — i *e > l (1) 

4ln •|Z| 

For sufficiently long series of measurements, the occurrence of the 
individual data points si of s1, s2, sT can be viewed as statistically 
independent. 

The probability P=P(s1, sT|p, Z) for an occurrence of all measured data 
1 5 points s1 , sT can accordingly be written as: 

T 

P{s w .. 9 s T \^Y)=Y\p{s t \^t)= 
t=i 

e z (2) 
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Here, the unknown variables, the mean value |u and the covariance Z, 
depend exclusively on a(brain) model that describes the measured data. 

The model assumes a linear statistical relationship between the individual BOLD 
signals: 

N 

5 s t = ^SyS j +8i for i=l, TV 
or 

s = Ss + e (3) 

where s designates the external influence on the individual BOLD signals, like a 
sensory input of sensory cells on the examined areas of the brain. 

10 The influence variables si and ej on different examined areas i and j can in 

this case be correlated throughout. 

The model parameters to be specified are consequently the coupling 
strengths Sj of the underlying coupling matrix S, the mean value jae of the external 
influence s and the covariance Xe of e. 

15 The mean value |a and the covariance S depend on these: 

Z = Z(« ff ) (4) 

In this known method of analysis the model parameters are then 
20 determined such that the probability P=P(s1, sT| given in (2) for the 
occurrence of the measured data is maximized. 

For this purpose, a known maximum likelihood estimation method 
(optimization) is applied such as described in T.W. Anderson, An Introduction to 
Multivariate Statistical Analysis, Chapter 3, John Wiley & Sons, Inc., New York, 
25 London, Sydney, 1994. 
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Using the relationships (4) in (2), an expression that is dependent on the 
coupling strengths Si, the mean value jae and the covariance Ee is obtained, which 
expression is maximized by the optimization. 

The optimization then leads to the desired coupling strengths S\ between 
5 the BOLD signals. 

These in turn enable detection of functional relationships between various 
areas of the brain in specific perception processes or motor tasks (functional 
connectivity). 

This known method of analysis exhibits the disadvantage that the 
10 measured fMRI signals are able to be modeled only insufficiently accurately or 
that the model is matchable only insufficiently accurately to the measured fMRI 
signals, and consequently the mode of operation or interaction of neuronal areas 
is only insufficiently mappable. This shortcoming could possibly lead to incorrect 
conclusions being made with regard to connective functionality. 

15 A software tool for fMRI analysis, an "fmri.pro", is known from Specification 

for "fmri.pro" software relating to quantitative fMRI analysis, obtainable on 
07.09.2001, under http://www.med.uni-uenchen.de/radin/html/ 
arbeitsgruppen/fmri/ccfmri.html. A device for implementing the fMRI technique is 
known from Specification of fMRI - device, obtainable on 07.09.2001, under 

20 http://www.unipublic.unizh.ch/campus/uninews/2001/0147/fmri.html. 

SUMMARY OF THE INVENTION 

An object of the present invention is to provide an improved method for the 
analysis of neuronal activities. The improved method of analysis is intended to be 
able to better explain measured fMRI signals and thus to better describe the mode 
25 of operation and interaction of neuronal areas than is the case with the above- 
described known method of analysis. 

This object is achieved in accordance with the invention by a method for 
the analysis of neuronal activities in neuronal areas using signals describing the 
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neuronal activities, wherein signals are determined with each signal describing the 
neuronal activity in one of the neuronal areas. 

A matchable coupling that is described by using matchable coupling 
variables that describe a statistical relationship between the matchably coupled 
5 signals forms the basis of all the signals, rather than just some of them. 

Probabilities for an occurrence of the signals are determined, whereby a 
statistical distribution forms the basis of the occurrence of the signals. 

All matchable coupling variables are determined by an optimization of the 
probabilities and hence matched. 

10 The neuronal activities are analyzed using the matchable coupling 

variables. 

The arrangement in accordance with the invention for the analysis of 
neuronal activities in neuronal areas by the use of signals describing the neuronal 
activities has units connected with one another which are configured such that 

- the signals can be determined, whereby each signal describes the neuronal 
activity in one of the neuronal areas, 

- a matchable coupling can form the basis of all, not just some of the signals, 
said matchable coupling being described by using matchable coupling 
variables which describe a statistical relationship between the matchably 
coupled signals, 

- probabilities for an occurrence of the signals can be determined, whereby a 
statistical distribution forms the basis of the occurrence of the signals, 

- all matchable coupling variables can be determined by an optimization of the 
probabilities and hence matched, 

- the neuronal activities can be analyzed using matchable coupling variables. 
It is essential to the invention that the matchable coupling, which is 

described by using the matchable coupling variables, forms the basis of all 
signals. By this means, absolutely all coupling variables are determined by the 
optimization of probabilities and hence matched. 
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The invention thus differs from the known method of analysis described 
above because in the known method matchable statistical coupling forms the 
basis of only some of the signals. 

Only these can be determined by means of the optimization of probabilities 
5 and hence matched. 

Thus, the known method of analysis presupposes a known, predetermined 
and stipulated neuronal structure. 

In contrast to this, no predetermined and fixed coupling structures are 
assumed in advance by the inventive approach. These emerge only within the 
10 context of the optimization. 

Through the interaction of an optimization method and a search method, 
i.e. the search for existing couplings and the determination of their optimum 
values, in the inventive approach both the coupling structure most probable on the 
basis of the signals and a coupling strength of the specific couplings are 
15 determined. 

A particularly advantageous aspect of the inventive approach is that this 
approach is independent of other methods and of possibly defective prior 
knowledge. No prior knowledge, or just outline prior knowledge, of coupling 
structures is needed in the case of the invention in order to analyze the neuronal 
20 activities. 

As a result of the flexibility achievable by the invention in the matching of 
couplings, neuronal structures can be determined with greater precision and in 
greater detail. 

The computer program with program code according to the invention is 
25 equipped to perform all steps in accordance with the method of analysis of the 
invention when the program is run on a computer. 
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The computer program product with program code stored on a machine- 
readable medium is equipped to perform all steps in accordance with the method 
of analysis of the invention when the program is run on a computer. 

The arrangement and the computer program with program code that are 
5 equipped to perform all the steps in accordance with the inventive method of 
program product are particularly suitable for performing the method of analysis 
according to the invention or one of its further developments explained below. 

The further embodiments described below relate both to the methods and 
to the arrangement. 

10 The invention and the further embodiments described below can be 

implemented both in software and in hardware, for example using a special 
electrical circuit (ASIC). 

Furthermore, implementation of the invention including the further 
embodiments described is possible by means of a computer-readable storage 
15 medium on which is stored the computer program with program code. 

The invention, including the further embodiments described below can also 
be implemented in a computer program product which has a storage medium on 
which is stored the computer program with program code which executes the 
invention. 

20 The statistical distribution that forms the basis of the occurrence of the 

signals can be a statistical distribution of a first or a higher order. The higher 
order can be achieved by using an Edgeworth expansion, as described in Samuel 
Kotz, Norman L Johnson (Editors-ln-Chief), Cornish-Fisher and Edgeworth 
Expansions, Ch. 4, pages 188-192, Encyclopedia of Statistical Sciences, Volume 

25 2, John Wiley & Sons, 1982 or a sum of normal distributions. 

In such a statistical distribution of a higher order not only do mean value 
and covariance have to be matched to comply with a quantity of data - as in the 
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case of such a statistical distribution of the first order - but further higher-order 
parameters like moments and cumulants also have to be matched. 

It should be noted that the possibilities mentioned for achieving a higher 
order without restricting the generality are only two selected statistical 
5 distributions. Other possibilities are known to persons skilled in the art. 

In addition, in the sum of normal distributions, the individual normal 
distributions and thus indirectly the neuronal activities can be weighted. 

Optimization can also be implemented by means of a maximum likelihood 
estimation, as described in T.W. Anderson, An Introduction to Multivariable 
10 Statistical Analysis, Chapter 3, John Wiley & Sons, Inc., New York, London, 
Sydney, 1994 method. 

In the optimization, a relationship between the linear statistical relationship 
and the statistical distribution can be taken into consideration as an auxiliary 
condition. 

15 It is also expedient, because the biological model of neuronal structures 

can be mapped in a more real way, for external influences on the signals to be 
taken into consideration in the linear statistical relationship. Such external 
influences may, for example, be sensory inputs of sensory cells on the examined 
areas. 

20 The determination of signals in the invention, for example of BOLD signals, 

can be carried out by measuring signals or by transferring and/or reading in 
existing signals. 

The invention is particularly suitable for use in an fMRI technique, which is 
considerably improved and more powerful as a result. 

25 In the context of such an fMRI use or fMRI examination, the neuronal areas 

are areas of the brain with corresponding nerve structures of patients to be 
examined and diagnosed. 
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In the fMRI examination using the inventive approach, BOLD signals are 
measured in various areas of the brain of a patient for defined perceptory or motor 
tasks carried out by the patient, which BOLD signals describe or represent the 
neuronal activities in the respective areas of the brain. These are evaluated or 
5 analyzed, whereby the signal coupling variables are determined. 

Using the analysis results, in particular the signal coupling variables, 
functional as well as physical dependencies between areas of the brain can be 
detected and determined. These can further be used for a diagnosis concerning a 
functional disorder in an area of the brain of a patient, for example by comparing 
1 0 "disordered" dependencies with those of healthy persons. 

DESCRIPTION OF THE DRAWINGS 

Figure 1 illustrates a functional magnetic resonance imaging apparatus for 
conducting an fMRI examination in accordance with the invention. 

Figure 2 is a flowchart for analyzing BOLD signals according to an 
1 5 exemplary embodiment of the invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The invention will be explained in the context of an embodiment for 
functional magnetic resonance imaging (fMRI). 

Fig. 1 shows an apparatus 100 for carrying out functional magnetic 
20 resonance tomography (fMRI). 

Basic principles of fMRI technology, which is a further development of the 
known modality of magnetic resonance tomography, are known from A. W. Toga 
and J. C. Maziotta (Ed.), "Brain Mapping: The Methods", Ch 9: M. S. Cohen: 
"Rapid MRI and Functional Applications", Academic Press 1996. 

25 The apparatus 100 has a closed tunnel 110 which is incorporated in a 

magnet 120 such that it generates a strong magnetic field in the tunnel 110. 
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The apparatus 100 also has a patient table 130 that can be introduced into 
the tunnel 1 10, on which table a patient is placed for an examination. 

In addition, the apparatus 100 has a control unit 131 that enables 
monitoring and control of the patient table 130 in the examination, for example a 
5 controlled introduction of the patient table 130 into the tunnel 120. 

The tomography apparatus 100 also has a measuring device 140 for the 
measurement of BOLD (Blood Oxygenation Level Dependent) signals, an 
associated evaluation device 141 for evaluating the measured BOLD signals, in 
this case a high-performance computer, as well as an operating or interface 
10 device 142 for operating personnel and a display device 143 for displaying an 
examination result. 

The components of the apparatus 100 are functionally connected to one 
another, for example via signal lines or data lines 150, via which the data and 
signals can be transferred. 

1 5 Based on the fMRI technique, the neuronal activity in areas of the brain of a 

patient can be measured and analyzed and a diagnosis can be made by means of 
the apparatus 100 shown in Fig. 1. 

To this end, the BOLD (Blood Oxygenation Level Dependent) signal in 
individual selected areas of the brain of the patient is measured by the measuring 
20 device 140, the BOLD signal being representation of the neuronal activity in the 
respective area. 

The result of such fMRI measurements shows the curve of the activity of 
the individual areas over a certain period of time, for example during cognitive 
sequences as a result of specific perception processes or motor tasks which have 
25 to be carried out by the patient during an examination. 

Functional disorders in the brain of the patient are thus inherently contained 
in the measured fMRI signals. 
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Using the evaluation device 141, which makes available or implements the 
inventive method of analysis, the fMRI measurements, i.e. the BOLD signals 
measured in individual areas of the brain, are analyzed. 

This inventive method of analysis represents an improvement over the 
5 known method of analysis described above. 

In the inventive method of analysis, brain activity is determined in the form 
of corresponding activation patterns in the examined areas in the brain and/or 
relationships between activation patterns in the examined areas and from that 
conclusions drawn directly of functional disorders in the brain and their causes. 

10 The inventive of analysis made available by the evaluation device 140 is 

based upon an extended and more flexible model of the brain, of the neuron 
structures in the brain and their behavior, in particular their interactions, on the 
basis of which the measured BOLD signal is analyzed and evaluated. 

The basic principles of the inventive method of analysis and the model are 
15 explained below. 

The results or the conclusions of an examination are shown on the display 
device 143 and can be further processed by means of the operating and interface 
device 142 in connection with the evaluation device 141. They also serve as a 
basis for a medical diagnosis of a patient to be examined and diagnosed. 

20 Basic principles of the inventive method of analysis as set forth in Fig. 2, 

steps 210 to 250. 

It should be noted that the inventive method of analysis is an improved 
further development of the known method of analysis described above. 
Consequently, - unless stated otherwise - the known and inventive method of 
25 analysis correspond in parts. Where matching parts are mentioned explicitly, they 
exhibit the labeling used previously hereinabove. 
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Using the inventive method of analysis 200 the fMRI measurements (210), 
i.e. the BOLD signals in examined areas of the brain of a patient, are analyzed 
(210 to 250) and/or compared with reference fMRI measurements. In this way, 
conclusions are drawn directly about functional disorders in the examined brain 
5 and their causes. 

The inventive method of analysis 200, which generates statistical 
characteristic quantities, such as statistical correlations between fMRI 
measurements in various areas of the brain, is based on an extended and more 
flexible mathematical model (220) of the brain based upon the known 
10 mathematical model according to A. W. Toga and J. C. Maziotta (Ed.), "Brain 
Mapping: The Methods", Ch 9: M. S. Cohen: "Rapid MRI and Functional 
Applications", Academic Press 1996. 

In this extended model (220) of the inventive method of analysis, the 
coupling matrix S is populated in all (matrix) positions by changeable coupling 
1 5 strengths Si. 

In the inventive method of analysis 200 all - because they are also 
changeable - coupling strengths Si are determined in such a way that statistical 
characteristic quantities which are determined from the fMRI measurements can 
best be explained (210 to 250). 

20 A data point s=s t represents the totality of all the BOLD signals s1 , sN of 

the individual n examined areas at a point in time t (or averaged over a time 
interval t) (t=[1;T]). 

The fMRI measurement comprises a large number of such data points s1, 
s2, sT for differing perception processes and/or motor tasks for which the 
25 corresponding BOLD signals were measured. 



In contrast to the known method of analysis, in which a multivariant 
standard distribution was assumed for the statistical distribution of the data points, 
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in the new method of analysis 200 a weighted sum of normal distributions is 
assumed for the statistical distribution (220). 

p(s|C 1 ,,c L ,/i 1 ,...// L ,5: 1 ,... s s L ) = 



5 equivalence of probabilities P=P(sl C1, CL, ji1, ^L, Z1, .... IL ) (230) (cf. 
(2)) for the occurrence of measured data points s1, s2, sT depend on more or 
different parameters than the mean value \x and the covariance X of the known 
method of analysis. 

In the inventive method of analysis 200 specific statistical variables, which 
10 can be calculated for the chosen statistical distribution, are now correlated with 
the model parameters, i.e. the coupling strengths Si, the mean value \jls of external 
influence e and the covariance Ze of s. 



ZL and all the moments and cumulants of the chosen distribution of a higher 
15 order. 

An inherent relationship between the statistical distribution parameters and 
the model parameters to be determined emerges from this, in this case taking into 
consideration the distribution (5) and the extended model based upon the model 
according to (3). 

20 ju = //(C l5 C L ,// 1 ,...,// Z/ ,E 1 ,...,i: L ) 

Z = I,(C l ,....,C L ,ju u ... 9 jU Ly Z u ...,I, L ) 



(5) 




In this case, the chosen statistical distribution and thus also the 



These include inter alia the means values ja1, j^L, the covariances Z1, 



(6) 
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In conformance with the known method of analysis, the optimum model 
parameters are now determined (240) in an analogous manner in the new method 
of analysis 200 using maximum likelihood estimation [1] through optimization or 
maximization of the probabilities (5). 

5 The basic principles of maximum likelihood estimation are described in 

T.W. Anderson, An Introduction to Multivariable Statistical Analysis, Chapter 3, 
John Wiley & Sons, Inc., New York, London, Sydney, 1994. 

The parameters to be taken into consideration for the optimization are the 
parameters of the chosen statistical distribution of a higher order, in this case the 
10 weighted sum of normal distributions, the desired model parameters and statistical 
variables, in this case the mean value m and the covariance Z from (6), via which 
the relationships between the model parameters and the statistical distribution (5) 
were produced. 

These relationships from (6) must be taken into consideration as auxiliary 
1 5 conditions in the optimization. 

The optimization then leads to the desired coupling strengths Sj that 
describe dependencies between the BOLD signals (250) and form the basis of the 
further evaluation and medical diagnosis (250). 

As an alternative, in place of the weighted sum of normal distributions, the 
20 distribution of data points can be described by an Edgeworth expansion. 

The basic principles of the Edgeworth expansion are described in Samuel 
Kotz, Norman L. Johnson (Editors-ln-Chief), Cornish-Fisher and Edgeworth 
Expansions, Ch. 4, pages 188-192, Encyclopedia of Statistical Sciences, Volume 
2, John Wiley & Sons, 1982. 

25 Although modifications and changes may be suggested by those skilled in 

the art, it is the invention of the inventors to embody within the patent warranted 
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heron all changes and modifications as reasonably and properly come within the 
scope of their contribution to the art. 

CH1\ 4226890.1 



